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SUMMARY 

A single elastic fiber embedded in an infinite elastic matrix is considered. An incident 
plane SH wave is assumed in the infinite matrix, and an expression is derived for the total 
energy in the fiber due to the incident SH wave. A nondimensional form of the fiber energy 
is plotted as a function of the nondimensional wavenumber of the SH wave. It is shown 
that the fiber energy attains maximum values at specific values of the wavenumber of the 
incident wave. The results obtained here are interpreted in the context of phenomena 
observed in acousto-ultrasonic experiments on fiber reinforced composite materials. 
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INTRODUCTION 


Quantitative ultrasonic nondestructive evaluation methods require analyses of wave 
propagation in the material or structure which is to be evaluated. In fiber reinforced com- 
posite materials, wave propagation characteristics depend on properties and arrangement 
of the constituent materials, as well as on the overall geometry and both the microscopic 
and macroscopic defect state of the composite. 

Analyses of wave propagation in fiber reinforced composites are reviewed in [1] and 
[2]. According to these reviews, wave propagation analyses may be divided into two broad 
categories. In the first category, the composite is treated as an equivalent homogeneous 
medium, and mechanical responses such as stress and strain are described in terms of 
averages over the constitutive elements. It has been stated that the homogeneous or 
equivalent continuum approach is valid “when the scale of the changes in stress level (rise 
distance, wavelength, etc.) is much larger than the sizes of the constituents of the composite 
(fiber or particle diameter, fiber spacing, ply spacing, etc.)” [1]. In the second category 
of wave propagation analyses, the composite is modeled as a heterogeneous mixture of 
constituents, and the dynamic interaction between the constituent materials is considered. 
The heterogeneous analyses are complex, but are necessary if the scale of changes in stress 
level is smaller than or of the same order as the size of the composite constituents, or if 
small-scale phenomena such as interface bond failure are of interest. 

Recent acousto-ultrasonic experiments by Kautz [3] show some interesting wave prop- 
agation behavior in graphite fiber reinforced epoxy composite panels. Kautz introduces a 
broadband signal into a composite panel, and observes that the response of the composite 
appears to be divided into a high frequency component traveling at a relatively fast wave 
speed and a low frequency component traveling at a relatively low wave speed. When the 
same broadband signal is introduced into a geometrically similar panel consisting only of 
the epoxy matrix, the response is observed to consist only of the low frequency component. 
Kautz hypothesizes that the composite panel acts as a heterogeneous medium, with wave 
propagation along the fibers at the higher wave speed and wave propagation through the 
matrix at the lower wave speed. 

In this report, an infinite elastic matrix containing a single cylindrical elastic fiber 
is considered. An incident plane SH wave is assumed to exist in the elastic matrix, and 
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an expression is derived for the total energy of the fiber due to the incident SH wave. 
It is shown that the total fiber energy attains maximum values at specific values of the 
wavelength of the incident wave. The purpose of this work is to investigate, in a simple 
setting, the issue of energy transfer between the matrix and the fibers in a composite 
material. Recommendations are given for future research which may give further insight 
into the experimentally observed phenomena in [3]. 

ANALYSIS 

A cylindrical fiber having circular cross section is shown in Fig. 1. Rectangular coordi- 
nates ( x,y,z ) and cylindrical coordinates (r, 6,z) are defined as shown in Fig. 1. Displace- 
ments u x > Uy and u z are defined in the x , y and z directions, respectively, and displacements 
u r , ue and u z are defined in the r, 0 and z directions, respectively. The fiber is assumed 
to have radius a and to be of infinite extent in the z direction. The fiber is assumed to be 
elastic, with mass density pf and Lame constants p/ and A/. It is assumed that the fiber 
is embedded in an infinite elastic matrix with mass density p and Lame constants p and 
A. Continuity of relevant stresses and displacements is assumed at the materials interface 
t — a. 

The infinite elastic matrix is assumed to contain an incident harmonic plane wave of 
the form 


i \x,y,z,t ) = 0 

(1) 

o 

II 

<<■— V 

(2) 

= Re lw 0 e 3 ( kx wt) j 

(3) 


where the superscript (*) denotes the incident wave, u>o is the amplitude of the incident 
wave, k is the wavenumber of the incident wave, u is radian frequency, t is time and 
j = y/—l. The wavenumber k and frequency u are related by 


< j = c s k 


(4) 


where 



(5) 
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Eqns. (l), (2) and (3) define a wave which propagates in the positive x direction with 
particle motion parallel to the fiber in the z direction. Such a wave, in this context, is 
denoted as an SH wave. 

Because of the cylindrical geometry of the fiber, it is convenient to write the incident 
SH wave in the cylindrical coordinates (r, 0, z ) as [4] 


u^(r,0,z,t) = 0 
u { e\r,0,z,t) =0 

t4 l) (r,0,z,f) = Re 


wq ^ e n j n J n (kr ) cos nde Jwt 


n— 0 


( 6 ) 

( 7 ) 

( 8 ) 


where J n is a Bessel function of the first kind and of order n, and where 




f 1, n=0 
\ 2, n=l,2,. . . 


( 9 ) 


It is shown in [4] that eqns. (1) through (3) and eqns. (6) through (8) are mathematically 
identical representations of the same incident wave. 

The interface conditions at r = a and the three-dimensional equations of elasticity 
which govern the infinite medium and the fiber are satisfied by the total displacement field 

[ 4 ] 


where 


u r (r,0,z,t) = 0 
Uff(r,0,z,t) = 0 

u z (r,0,z,t) = u^(r,0,z,t) + u ( z g) (r,0,z,t ) 


r > a 


u r (r,0,z,t ) = 0 
u g (r,0,z,t) = 0 
u z (r,0,z,t) = ui f) (r,0,z,t) 


r < a 


u^( r ,0,z,t) = Re < w 0 ^ A n H^\kr) cos n0e ,wt 


n=0 


u i^( r , 0,2,0 = Re 


| -w 0 jr C n J n 

l n=0 


(k /r) cos n0e ]wt 


( 10 ) 

( 11 ) 

( 12 ) 

(13) 

(14) 

(15) 

(16) 
(17) 
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A n = (-j n e n /A)[nfkfJ n {ka)J' n (kfa) - ^kJ' n {ka)J n {k f a)\ 
Cn = {-j n tn/A)nk J' n (ka) H ^ (ka) - {ka) J n {ka ) j 

A = p f k f HW{ka)J’ n (k f a) - /ikl!™' {ka)J n (k f a) 


kf = -^~k 
c *f 


The quantities A n and C n are dimensionless complex coefficients. is a Hankel function 
of the first kind and of order n, and the prime denotes the derivative of a function with 
respect to its argument. 

Eqn. (16) represents a scattered wave which propagates outward from the fiber. As 
r — > oo, the amplitude of the scattered wave approaches zero, so that far from the fiber 
the total displacement field approaches the displacement field of the incident SH wave. 

Eqn. (17) represents a refracted or transmitted wave in the fiber. Eqn. (17) is a 
standing wave with time-harmonic particle motion in the z direction. Over any fiber cross 
section z — constant, the fiber displacement in the z direction varies as a function of r and 
0 according to eqn. (17). 

The nonzero fiber strains associated with the displacements given by eqns. (13) through 
(15) are 

1 du * /) 


1 dui P 


r 36 

duP 


Irz = 


and the nonzero stresses in the fiber are 


OBz = 


r 80 

3u i'> 


= w-5T 

The elastic strain energy U j per unit length of fiber in the z direction is 

ra /*2ir i 


/ a rzir j 

J - [Orzlrz + oezlez) r dOdr 
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The kinetic energy K / per unit length of fiber in the z direction is 


ra r2ir -■ 2 

Kf = Jo Jo 2 Pf \ z J rd6dr 

The total energy Ef per unit length of fiber is, therefore, 


( 28 ) 


Ef = U f + Kf 


(29) 


The energy Ef given by eqn. (29) is a function of time. In order to obtain a single 
numerical measure of energy, it is customary to average the energy Ef over one time period 
of the harmonic oscillation. The time-averaged energy < E > / is defined by 


< E >f 




dt 


where 


r = 


2n 


U) 


(30) 


(31) 


When eqn. (17) and eqns. (23) through (29) are substituted into eqn. (30), a long 
algebraic calculation gives the following result for the time-averaged fiber energy: 


< 


+J*(kfa) - J n -i(kfa)Jn+i(k f a) J ( 32 ) 

r n (*/“)]} 


+nJi 


where 


= / 1/2, n = 0 

\l, n = 1, 2, . . . 


(33) 


and where the asterisk denotes the complex conjugate. The derivation of eqn. (32) requires 
the complex variable identity 

(34) 


= ~(z + z*) , 


the trigonometric identies [5] 


L 


2k 


cos nB cos mO dO = —8 mn 

0 Pn 


(35) 
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f 2w f 0 

/ sin n0 sin mO — < 

Jo l* 

and the Bessel function identities [6] 


for n=0 or m=0 
ir6 mn otherwise 


(36) 


Q r ^n(^f r ) — k / d n+ i{kjr) + ^ J n (kjr) 


(37) 

(38) 

(39) 


rJ l( k f T ) dr = ^ - Jn-i(k/a)J„+i(k/a)} 

J q Jn{k f r)J n +i{k f r)dr = i- | i [l - J 0 2 (fc/o)] - ^ Jj{k f a) 

J q \jl{ k f r ) dr = ^ |i + J oi k f a ) + J ni k f a ) J *(*/ a ) J > n >° ( 40 ) 

The quantity d mrt in eqns. (35) and (36) is the Kronecker delta. 

Eqn. (32) is an expression for the total fiber energy per fiber length in the z direction. 
In order to nondimensionalize the fiber energy < E > f, consider an infinite uniform elastic 
matrix containing only the incident SH wave given by eqns. (6) through (8), and consider 
the time-averaged energy < E > m contained, per unit length in the z direction, in the 
region 0 < r < a. By a calculation similar to the calculation leading to eqn. (32), the 
time-averaged energy < E > m is 


< E > m = \nwl f^{D* n D n j- ^{j 2 +1 (M - J n (ka)J n+ 2 (ka) 


(ka)‘ 


d" Jn—l{kd)Jn+l{kO')'^'^~TlJ n ^ka) ^ 


where 


Dn = Cni n 


(41) 

(42) 


A nondimensional fiber energy < e >/ is now defined by 


< e >/= 


< E >/ 
<E> m 


(43) 


The dimensionless quantity < e >f gives a measure of the fiber energy relative to the 
energy the matrix would contain if there were no fiber present. 
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By introducing the nondimensional parameters 


r P = 


r M = 


£/ 

_ HL 


(44) 

(45) 


eqns. (43) and (19) can be written as 


<e >/= 


<£'>/ 


< E >, 


t-^co j_ (nffca) 2 f Jn+i(~lka) J n (-)ka)J n+ 2 (^ka) 1 , jl \ 

2.n=0 2 | + J2 {lka) _ J n _ , (^d) J n+1 (^ka) / + nJ ^ ka ) 

— T u r s _ o , - _ , . v _ , . v x -i 


Y,™=0 D n D n L 


(fc°) 2 / ^n+i(^ a ) Jn{ka)J n + 2{ka) 1 r2('t /7 'i 

2 \ +Jl{ka) - ./^(/ca) J n+ 1 (A;a) / + nJ n W 


(46) 


C n = 


{-j n e n ) [j^{ka)H { n 1 \ka) - H { n iy {ka)J n {ka) 
r^H^\ka)J^ka) - H^ iy (ka)J n (')ka) 


respectively, where 


1 


-te) 


1/2 


(47) 


(48) 


From eqns. (46), (47) and (48) it can be seen that the nondimensional fiber energy < e >/ 
depends only on the nondimensional parameters ka, r p and r M . The parameter ka is a 
nondimensional wavenumber of the incident wave, and can be written in terms of the 
incident wavelength A,- as 


ka = 


27 ra 

"*7 


or in terms of frequency as 


, ua 
ka= — 
c 

The nondimensional fiber energy < e > / is plotted as a function of the nondimensional 
wavenumber ka in Fig. 2 for the numerical values r p = 2.0 and r M = 25.0. These numerical 
values correspond to an E-glass fiber in a PMR (polyimide monomeric reactant) matrix 
[7]. In Fig. 3, the nondimensional fiber energy < e >/ is plotted as a function of ka for the 
numerical values r p = 1.5 and r M = 7.5. These numerical values correspond to a Celion 
6000 fiber in a PMR matrix [7]. The numerical values used to generate Fig. 3 are estimates 
of the material properties of the graphite/epoxy composite used in [3]. The infinite series 
in eqn. (46) are computed by truncating when the absolute value of the individual terms 
becomes less than 10~ 12 times the absolute value of the accumulated sum. 
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DISCUSSION 


In both Fig. 2 and Fig. 3, the nondimensional fiber energy < e > / is approximately 
equal to one when ka is close to zero, or, equivalently, when the wavelength of the incident 
wave is long compared to the fiber diameter. This means that for long wavelengths, the 
energy in the fiber is approximately the same as the energy which would exist in a uniform 
matrix. In other words, the presence of the fiber does not significantly affect the energy 
distribution for long wavelengths. As ka becomes larger (or, equivalently, as the wavelength 
becomes smaller) the energy in the scattered wave becomes larger, and the energy in the 
fiber becomes smaller than the energy which would exist in the region 0 < r < a in the 
uniform matrix. 

The most interesting features of Fig. 2 and Fig. 3 are the local maxima in the values 
of < e >f. For an incident wave containing a broad band of wavenumber, the fiber energy 
will be greatest at frequencies corresponding to the values of ka where the local maxima 
occur. If the fiber response could be measured, the local maxima would appear in a 
frequency spectrum of the measured response. Thus it is sensible, at least in this simple 
case, to suppose that certain frequency components in a measured response may contain 
information about the fiber. 

The incident SH wave assumed here is independent of the spatial coordinate z and so, 
therefore, is the displacement field in the fiber. Therefore, it is impossible, in this analysis, 
to discuss the issue of wave propagation along the fiber in the z direction. The analysis of 
an incident wave which has a z component may allow consideration of wave propagation 
along the fiber. Equations which describe the scattering and refraction of such an obliquely 
incident wave are given in [4]. 

CONCLUSIONS AND RECOMMENDATIONS 

The purpose of this report is to investigate the transfer of energy to a single fiber in 
an infinite elastic matrix containing an incident SH wave. An expression is derived for the 
energy in the fiber, and it is shown that the fiber energy attains local maxima at specific 
values of the wavenumber of the incident wave. 

Future research should consider the problems of incident P and SV waves. As men- 
tioned above, an obliquely incident wave may allow analysis of wave propagation along 
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the fiber. Also, properties of the scattered wave given by eqn. (16) should be investigated, 
since the scattered wave may be available for experimental measurement. It would be in- 
teresting to see if the energy flux of the scattered wave has wavelength dependence similar 
to the wavelength dependence of the energy of the fiber. 
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Fig. 1 Elastic fiber embedded in infinite elastic matrix containing incident SH wave. 
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Fig. 2 Nondimensional fiber energy, < t > f as a function of nondimensional wavenumber ka 
for r p = 2.0, r M = 25.0. 
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Fig. 3 Nondimensional fiber energy, < e >y as a function of nondimensional wavenumber ka 
for r p = 1.5, and r M = 7.5. 
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